Method and devices for synchronization

ABSTRACT

This invention relates to methods and devices for time and frequency synchronization. The invention has particular application where time and frequency synchronization over packet networks using, for example, the IEEE 1588 Precision Time Protocol (PTP) is being carried out. The primary challenge in clock distribution over packet networks is the variable transit delays experienced by timing packets, packet delay variations (PDVs). Embodiments of the invention provide a method for time offset alignment with PDV compensation where a synchronized frequency signal is available at a slave device via Synchronous Ethernet and is used to determine the compensation parameters for the PDV.

FIELD OF THE INVENTION

The present invention relates to methods and devices for synchronization. It is particularly, but not exclusively, concerned with time synchronization over packet networks using, for example, the IEEE 1588 Precision Time Protocol (PTP) and synchronous Ethernet.

BACKGROUND OF THE INVENTION

There are many applications of high precision time synchronization over packet networks based on packet technologies like IP and Ethernet. Synchronization is a key requirement in telecommunication networks and the telecommunication industry has become a driving force in the development and evolution of synchronization solutions and standards. For example, wireless technologies like GSM, UMTS (both frequency division duplexing (FDD) and time division duplexing (TDD) versions of the UMTS technologies), and CDMA2000 require frequency accuracy of 0.05 ppm (parts per million) at the air interface.

In addition, CDMA2000 requires time synchronization at the ±3 μs level (±10 μs worst case) and UMTS TDD modes (i.e., TD-CDMA and TD-SCDMA) require accuracy off 1.25 μs between base stations. The accurate reference clock is typically derived from Time Division Multiplexing (TDM) interfaces (for GSM, UMTS-FDD (i.e., W-CDMA) and other FDD based technologies) or from expensive GPS receivers located at the base station (for TDD based technologies). In reality, TDD base stations today use GPS and wireless operators have long sought alternative synchronization sources.

Without synchronization traceable to a highly accurate reference clock in the wireless network, local interference between channel frequencies, as well as mutual interference with neighboring base stations will occur. This could ultimately cause calls to be dropped and a degradation of the overall user experience. The high level of synchronization provided by IEEE 1588 PTP could also be useful in time-of-arrival (TOA) and time-difference-of-arrival (TDOA) based localization solutions where one is interested in locating mobile nodes connected to a wireless network.

IEEE 1588 PTP has been designed as an improvement to current methods of synchronization within a distributed network of devices connected to a packet network. It is designed for systems requiring very high accuracies beyond those attainable using Network Time Protocol (NTP). Unlike NTP, IEEE 1588 uses hardware timestamping to deliver timing accuracy well under a microsecond (nanosecond level accuracies). It is also designed for applications that cannot bear the cost of a GPS receiver at each node, or for which GPS signals are inaccessible.

IEEE 1588 PTP can be used for distributing timing information (frequency and time) from a master to one or more slaves. The clock recovery mechanism at a slave will normally have to correct for both frequency and time deviations from the master clock during clock recovery.

Synchronous Ethernet (Sync-E), defined in a suite of ITU-T Recommendations, describes a method for frequency (only) distribution over the physical layer of Ethernet (IEEE 802.3). Since time synchronization requires both frequency and time offset alignment between master and slave clocks, some solutions suggest the combine use of Sync-E and IEEE 1588 PTP where Sync-E provides frequency alignment while IEEE 1588 PTP provides time offset correction.

In this setup of hybrid Sync-E and IEEE 1588 PTP, Sync-E provides perfectly a frequency synchronization signal using Ethernet physical layer frequency transfer, while IEEE 1588 PTP provides time synchronization over the packet network (above the physical layer). With frequency synchronization provided, the only issue here is determining the time offset between slave and master clocks. The time offset should be fixed (a constant) since slave's frequency is aligned to the master's. However, even for time offset determination, the problem here is that IEEE 1588 PTP messages are still exposed to the PDVs inherent in the packet network even though frequency synchronization is available through Sync-E. This is particularly the case if PTP messages are sent end-to-end from master to slave with the intermediate packet network providing no assistance whatsoever (like hop-by-hop Boundary Clocks or Transparent Clocks) to mitigate and compensate for PDV.

PDV therefore still makes time synchronization a challenge even though the frequency synchronization aspect is solved.

A complete IEEE 1588-based solution includes clock recovery (control) algorithms, filters, and PTP-clock based on hardware timer and direct timer access. The IEEE 1588 standards define a wide range of synchronization capabilities except the clock recovery mechanisms (servo algorithm, phase-locked loop (PLL), timers, etc.) to be used at the receiver (slave) to synchronize its local clock to the master. The last pieces are vendor and proprietary solutions and are often product differentiators.

Effects of Packet Delay Variations on Clock Synchronization

PDV is the main factor affecting the accuracy and stability of slave clocks when using packet timing protocols such as IEEE 1588 PTP. Packet network devices such as switches and routers introduce a variable delay to packets (PDV) that inhibits accurate path delay measurements and clock synchronization. Even for techniques that do not require path delay measurements for clock (e.g., frequency) synchronization, the PDV is a direct contributor to the noise in the recovered clock. The PDV inherent in packet networks is a primary source of clock noise. The higher the clock noise, the poorer the clock quality rendering the recovered clock sometimes unusable for end system applications when the noise exceeds application defined thresholds. (The term clock noise here refers to all impairments to the timing information recovered at the slave including jitter, wander, and other imperfections in the recovered clock.)

For instance, for time synchronization, the variation in delay from packet to packet through the network induces noise in the slave's perception of the time at the master. Constant delay (assuming frequency synchronized clocks) would cause a fixed offset, however variable delay causes a varying estimate of the offset. The performance of the slave is affected by both the magnitude of this variation, and how effective the slave's filter is at removing this noise.

Synchronous Ethernet

Clock synchronization methods fall under two broad categories, namely: Layer 1 methods, and Layer 2 and higher methods (Layer 2+) of the OSI model of communication. A schematic of the Layer 1 methods is shown in FIG. 1. A schematic of the Layer 2+ methods is shown in FIG. 2.

The Layer 1 (e.g., Sync-E, PDH, SONET/SDH), and Layer 2+ methods (e.g., NTP, IEEE 1588 PTP) differ mainly in the method used to inject and transport synchronization over the network 2 from a reference clock 1. Layer 1 methods achieve synchronization via the physical layer 4a, 4b bit stream while the Layer 2+ methods do so using higher levels 5a, 5b above the physical layer, for example, via differential clocking mechanisms, processing transmitted clock samples (timestamps) encoded within the packets from a transmitter to a receiver, or processing packet arrival patterns at the receiver to generate timing signals. Layer 1 methods generally support frequency transfer only while Layer 2+ methods support both frequency and time transfer.

ITU-T Recommendations G.8261, G.8262, G.8264 (which define Sync-E) indicate that synchronization can be achieved over Ethernet links by synchronizing the bit clock of the physical layer as is currently done on SONET/SDH links. The clock distribution aspect of Sync-E is essentially native Ethernet equipped with point-to-point distribution of timing signals. By enabling point-to-point distribution of timing signals from an accurate timing reference source, Sync-E aims to bring carrier grade telecom quality clocks to packet networks. Sync-E, unlike IEEE 1588, cannot deliver absolute time (i.e., time-of-day, wall-clock); Sync-E provides frequency distribution only.

The most significant aspect of native Ethernet communications is that the transmitter and receiver clocks are independent (free-running with frequency accuracy not exceeding ±100 parts per million (ppm)) as shown in FIG. 3 and are not synchronized as in TDM networks. In fact, there need be no strict timing relationship between successive frames as individual frames are separated by an idle period or at a minimum an inter-packet gap (IPG). The Ethernet transmitter delimits each packet by a start sequence (“preamble”) and a stop sequence (“start-of-idle” or “end-of-stream delimiter”). At the receiver, a clock of the same nominal frequency is recovered and used to clock-in the data to a receive shift register. Only data that are bounded by the correct start and stop bit patterns are accepted.

Sync-E implements timing similar to what is seen in traditional TDM networks (e.g., PDH, SONET/SDH): the physical layer of an Ethernet link is used to deliver timing from one end of the link to the other. Timing transfer is achieved by enslaving the downstream link's Ethernet physical layer clock to the upstream clock. FIG. 4 shows the transmission of frequency from a master 10 to a slave 20. Synchronization is achieved through the use of a digital phase-locked loop (PLL) 6 a, 6 b in each of the master and slave. The PLL cleans (improves) the incoming clock signal in the slave 20 to filter out jitter and wander and in turn generates the necessary output clock(s), telecom and/or Ethernet clocks, for the receiving interfaces and devices.

FIG. 5 shows the way in which timing transfer is achieved across a sequence of Ethernet nodes 30 a-c (although 3 nodes are shown in FIG. 5, it will be appreciated that this configuration can be extended to cover significantly larger networks). Each node has its own crystal oscillator 31 which drives a local clock. The local clock is synchronized to the reference frequency 1 by timing transfer over the physical layer 4 components and the clock signals recovered using the PLL 6 a-c as described above. Effectively each node acts as a slave to its neighbour node in the “upstream” direction of the node 30 a having the frequency reference 1 and as a master to its neighbour node(s) in the “downstream” direction.

The obvious advantage of this approach is that, by relying on the physical layer only, clock synchronization quality is not influenced by impairments introduced by the upper layers, namely, packet delay variation (PDV), packet losses, and out-of-order arrival of packets at their destinations.

The 10 Mb/s version of Ethernet cannot support Sync-E. In fact, the 10 Mb/s Ethernet (10BASE-T) is not even capable of synchronization signal transmission over the physical layer interface because a 10BASE-T transmitter stops sending pulses during idle periods. A 10BASE-T transmitter simply sends a single pulse, keepalive pulse, every 16 ms to notify its presence to the receiving end. Of course, such infrequent pulses are not sufficient for clock recovery at the receiver. Idle periods in faster Ethernet flavors (100 Mb/s, 1 Gb/s and 10 Gb/s) are continuously filed with pulse transitions, allowing continuous high-quality clock recovery at the receiver and are therefore good candidates for Sync-E.

IEEE 1588 PTP

IEEE 1588 is now the industry accepted packet-based method/standard for synchronization of clocks in distributed systems with accuracies in the nanosecond levels. PTP can be used for both time and frequency synchronization, unlike Sync-E which allows for only frequency synchronization. PTP is a message based protocol that can be implemented across packet based networks including, but not limited to, Ethernet. Its underlying principle is a master/slave concept based on the regular exchange of synchronization messages. IEEE 1588 synchronizes all clocks within a network by allowing clocks to adjust their frequency/time to the highest quality clock (the GrandMaster clock).

IEEE 1588 PTP allows for two different types of timestamping methods, either one-step or two-step. One-step clocks update time information within event messages (Sync and Delay_Req) on-the-fly, while two-step clocks convey the precise timestamps of packets in general messages (Follow_Up and Delay_Resp).

A Sync message is transmitted by a master to its slaves and either contains the exact time of its transmission or is followed by a Follow_Up message containing this time. In a two-step ordinary or boundary clock, the Follow_Up message communicates the value of the departure timestamp for a particular Sync message. A Delay_Req message is a request for the receiving node to return the time at which the Delay_Req message was received, using a Delay_Resp message.

The message exchange pattern for the two-step clock can be explained as follows and is illustrated in FIG. 6. The master 101 sends a Sync message to the slave 103 over the packet network 102 and notes the time T₁ at which it was sent. The slave 103 receives the Sync message and notes the time of reception T₂. The master 101 conveys to the slave the timestamp T₁ by one of two ways: 1) Embedding the timestamp T₁ in the Sync message (a one-step clock—not illustrated). This requires some sort of hardware processing (i.e., hardware timestamping) for highest accuracy and precision. 2) Embedding the timestamp T₁ in a Follow_Up message (the two-step clock as shown in FIG. 6). Next, the slave 103 sends a Delay_Req message to the master 101 and notes the time T₃ at which it was sent. The master 101 receives the Delay_Req message and notes the time of reception T₄. The master 101 conveys to the slave 103 the timestamp T₄ by embedding it in a Delay_Resp message.

At the end of this PTP message exchange, the slave 103 possesses all four timestamps {T₁, T₂, T₃, T₄}. These timestamps may be used to compute the offset of the slave's clock 105 with respect to the clock reference 104 in the master 101 and the mean propagation time of messages between the two clocks. The computation of offset and propagation time assumes that the master-to-slave and slave-to-master propagation times are equal—symmetrical communication path.

The present invention aims to provide a method for time offset alignment with IEEE 1588 PTP with PDV compensation.

SUMMARY OF THE INVENTION

An exemplary embodiment of the invention provides a method of synchronising the frequency and time of a slave clock in a slave device to a master clock in a master device, wherein the master device and the slave device are connected by an Ethernet network, the method including the steps of: synchronising the frequency of the slave clock to that of the master clock using Synchronous Ethernet communications between the master device and the slave device; receiving in the slave device first and second timing messages sent from the master device and respective first and second timestamps being the time of sending of the first and second messages according to the master clock, and recording the time of receipt of the first and second messages according to the slave clock; calculating the forward variable delay encountered by the first and second messages in travelling from the master device to the slave device over the network by comparing the difference between the interval between the time of sending of said messages and the interval between the time of receipt of said messages; receiving in the master device third and fourth timing messages sent from the slave device and recording the time of receipt of the third and fourth messages according to the master clock; recording in the slave device the time of sending of the third and fourth messages according to the slave clock, sending the recorded time of receipt of the third and fourth messages to the slave device; calculating the reverse variable delay encountered by the third and fourth messages in travelling from the slave device to the master device over the network by comparing the difference between the interval between the time of sending said messages and the interval between the time of receipt of said messages; calculating, from the forward and reverse variable delays, a measure of the noise associated with both the receipt of the first and second timing messages and receipt of the third and fourth timing messages; estimating the variance of the measurement noise; and adjusting the actual time of the slave clock to synchronize it with the time of the master using the first and third timestamps, the recorded times of receipt of the first and third messages and the estimated variance.

A further exemplary embodiment of the invention provides a slave device connected to a master device having a master clock over an Ethernet network, wherein the slave device includes: a slave clock; an Ethernet adapter; and a processor, wherein: the Ethernet adapter is arranged to receive, at the physical layer, timing signals from the master device via Synchronous Ethernet and the processor synchronises the frequency of the slave clock to the frequency of the master clock using those timing signals; and the Ethernet adapter is arranged to receive, through a layer above the physical layer, first and second timing messages sent from the master device and respective first and second timestamps being the time of sending of the first and second messages according to the master clock, and to record the time of receipt of the first and second messages according to the slave clock; and further wherein the processor is arranged to: calculate the forward variable delay encountered by the first and second messages in travelling from the master device to the slave device over the network by comparing the difference between the interval between the time of sending of said messages and the interval between the time of receipt of said messages; send third and fourth timing messages to the master device for the master device to record the time of receipt of the third and fourth messages according to the master clock; record the time of sending of the third and fourth messages according to the slave clock; receive, from the master device, the recorded time of receipt of the third and fourth messages; calculate the reverse variable delay encountered by the third and fourth messages in travelling from the slave device to the master device over the network by comparing the difference between the interval between the time of sending said messages and the interval between the time of receipt of said messages; calculate, from the forward and reverse variable delays, a measure of the noise associated with both the receipt of the first and second timing messages and receipt of the third and fourth timing messages; estimate the variance of the measurement noise; and adjust the actual time of the slave clock to synchronize it with the time of the master using the first and third timestamps, the recorded times of receipt of the first and third messages and the estimated variance.

A further exemplary embodiment of the present invention provides an apparatus for synchronising the frequency and time of a slave clock in a slave device to a master clock in a master device, the apparatus including: the master device having the master clock and a first Ethernet adapter; the slave device having the slave clock, a second Ethernet adapter and a processor; and an Ethernet network connecting the master device and the slave device, wherein: the first Ethernet adapter is arranged to send, at the physical layer, timing signals based on said master clock using the Synchronous Ethernet protocol; the second Ethernet adapter is arranged to receive, at the physical layer, said timing signals from the master device and the processor synchronises the frequency of the slave clock to the frequency of the master clock using those timing signals; the master device is arranged to send, through the first Ethernet adapter, using a layer above the physical layer, first and second timing messages and respective first and second timestamps being the time of sending of the first and second messages according to the master clock; the second Ethernet adapter is arranged to receive, through a layer above the physical layer, said first and second timing messages and said first and second timestamps, and to record the time of receipt of the first and second messages according to the slave clock; the slave device is arranged to send, through the second Ethernet adapter, using a layer above the physical later, third and fourth timing messages, and to record the time of sending of the third and fourth messages according to the slave clock; the first Ethernet adapter is arranged to receive, through a layer above the physical layer, said third and fourth timing messages and sent from the slave device and to record the time of receipt of the third and fourth messages according to the master clock and send the recorded time of receipt of the third and fourth messages to the slave device; and further wherein the processor in the slave device is arranged to: calculate the forward variable delay encountered by the first and second messages in travelling from the master device to the slave device over the network by comparing the difference between the interval between the time of sending of said messages and the interval between the time of receipt of said messages; calculate the reverse variable delay encountered by the third and fourth messages in travelling from the slave device to the master device over the network by comparing the difference between the interval between the time of sending said messages and the interval between the time of receipt of said messages; calculate, from the forward and reverse variable delays, a measure of the noise associated with both the receipt of the first and second timing messages and receipt of the third and fourth timing messages; estimate the variance of the measurement noise; and adjust the actual time of the slave clock to synchronize it with the time of the master using the first and third timestamps, the recorded times of receipt of the first and third messages and the estimated variance.

BRIEF DESCRIPTION OF THE DRAWINGS

Embodiments of the invention will now be described by way of example with reference to the accompanying drawings in which:

FIG. 1 shows a schematic clock distribution method in which clock synchronisation information is transmitted over physical layer Ethernet and has already been described;

FIG. 2 shows a schematic clock distribution method in which clock synchronisation information is transmitted over upper layer Ethernet and has already been described;

FIG. 3 shows the distribution of timing information in native Ethernet and has already been described;

FIG. 4 shows the architecture of a synchronous Ethernet link between a master and a slave and has already been described;

FIG. 5 shows the distribution of timing information using synchronous Ethernet and has already been described;

FIG. 6 shows the message flow in a two-step clock under IEEE 1588 PTP and has already been described;

FIG. 7 shows, in schematic form, the distribution of timing information according to an embodiment of the present invention;

FIG. 8 shows simple models of the effect of clock skew on clocks with and without offset;

FIG. 9 shows the distribution and processing of timing information in a slave device according to an embodiment of the present invention;

FIG. 10 shows the effect of packet delay variations in the arrival of PTP messages at a slave device;

FIG. 11 shows how the packet delay variations can be measured using a synchronous Ethernet clock according to an embodiment of the present invention;

FIG. 12 shows the processing of timing information in a slave device according to an embodiment of the present invention; and

FIG. 13 shows how the slave clock can be adjusted to synchronize it with the master clock.

DETAILED DESCRIPTION

Accordingly, at their broadest, methods of the present invention provide for synchronising the time and frequency of a clock in a slave device to the time and frequency of a master clock in a master device by using both Synchronous Ethernet communications and additional timing messages.

A first aspect of the present invention provides a method of synchronising the frequency and time of a slave clock in a slave device to a master clock in a master device, wherein the master device and the slave device are connected by an Ethernet network, the method including the steps of: synchronising the frequency of the slave clock to that of the master clock using Synchronous Ethernet communications between the master device and the slave device; receiving in the slave device first and second timing messages sent from the master device and respective first and second timestamps being the time of sending of the first and second messages according to the master clock, and recording the time of receipt of the first and second messages according to the slave clock; calculating the forward variable delay encountered by the first and second messages in travelling from the master device to the slave device over the network by comparing the difference between the interval between the time of sending of said messages and the interval between the time of receipt of said messages; receiving in the master device third and fourth timing messages sent from the slave device and recording the time of receipt of the third and fourth messages according to the master clock; recording in the slave device the time of sending of the third and fourth messages according to the slave clock; sending the recorded time of receipt of the third and fourth messages to the slave device; calculating the reverse variable delay encountered by the third and fourth messages in travelling from the slave device to the master device over the network by comparing the difference between the interval between the time of sending said messages and the interval between the time of receipt of said messages; calculating, from the forward and reverse variable delays, a measure of the noise (“the measurement noise”) associated with both the receipt of the first and second timing messages and receipt of the third and fourth timing messages; estimating the variance of the measurement noise; and adjusting the actual time of the slave clock to synchronize it with the time of the master using the first and third timestamps, the recorded times of receipt of the first and third messages and the estimated variance.

The measure of the noise calculated may be referred to herein as “the measurement noise” but need not be an exact measure of the measurement noise in the system. The estimated variance may be an estimate of the variance of the measure of the noise obtained, which is used as an estimate of the variance of the actual measurement noise, or an estimate of the variance of the actual measurement noise which is produced using the measure of the noise.

The method of this aspect can provide a new approach for time synchronization where the frequency synchronization features of Synchronous Ethernet (Sync-E) are used by a slave device to extract certain packet transport measurement parameters which can then be used as an input to a time offset computation mechanism used for clock alignment at the slave.

Time synchronization requires both frequency and time offset alignment between a master and slave clock. In the method of the present aspect, frequency synchronization is provided via Sync-E, so the remaining computation is for the slave clock to determine its time offset from the master clock.

The time offset should be fixed (a constant) since slave's frequency is aligned to the master's. However, even for time offset determination, the problem is that timing messages sent end-to-end over the packet network are still exposed to the packet delay variations (PDVs) inherent in the packet network even though frequency synchronization is available through Sync-E. The PDV inherent in packet networks is a primary source of clock noise. The higher the clock noise, the poorer the clock quality rendering the recovered clock sometimes unusable for end system applications when the noise exceeds application defined thresholds.

Preferably the timing messages are messages according to the IEEE 1588 Precision Time Protocol, and more preferably the first and second timing messages are Sync messages and the third and fourth messages are Delay_Req messages.

The IEEE 1588 Precision Time Protocol (PTP) is a method for distributing high precision timing information (both frequency and time) from a master to one or more slaves over a packet network. In a conventional IEEE 1588 PTP arrangement, the slave then uses a clock recovery mechanism to recover the master clock from this timing information. The primary challenge in clock distribution over packet networks is dealing with the variable transit delays experienced by timing packets, that is, the packet delay variations (PDVs). The clock recovery mechanism at a slave will normally have to correct for both frequency and time deviations from the master clock during clock recovery.

Synchronous Ethernet (Sync-E), defined in a suite of ITU-T Recommendations, provides a method for frequency (only) distribution over the physical layer of Ethernet (IEEE 802.3). Since time synchronization requires both frequency and time offset alignment between master and slave clocks, the method of the present aspect may combine use of Sync-E and IEEE 1588 PTP where Sync-E provides frequency alignment while IEEE 1588 PTP provides time offset correction.

In this arrangement, Sync-E provides a perfect frequency synchronization signal using Ethernet physical layer frequency transfer, while IEEE 1588 PTP provides time synchronization over the packet network (above the physical layer).

However, the IEEE 1588 PTP messages (or other timing messages) are still exposed to the PDVs inherent in the packet network even though frequency synchronization is available through Sync-E. These PDVs still make time synchronization a challenge even though the frequency synchronization aspect is solved. In preferred implementations of the method of this aspect, time offset alignment with IEEE 1588 PTP with PDV compensation is provided where the synchronized frequency signal available at the slave via Sync-E is used to determine the PDV compensation parameters.

By combining features of Sync-E and IEEE 1588 PTP, the method of this aspect may allow the slave to accurately measure the actual PDV of PTP messages arriving at the slave. The slave can then correctly adjust its clock while compensating for the actual delay variations. This may allow the slaves to remove the negative effects that these delay variations cause.

In one implementation of the present aspect, the step of adjusting the actual time of the slave clock uses a linear Kalman filter to update a state equation representing the offset and skew of the slave clock compared to the master clock using the estimated variance of the measurement noise following receipt of each timing message in the slave device.

Practical implementation of the Kalman filter typically requires getting a good estimate of the process noise and measurement noise covariances. In the method of the present aspect, the frequency synchronization features of Sync-E can be used to extract accurately the PDV parameters from the PTP message flows to the slave. The PDV parameters can then be used to determine the measurement noise covariances for the Kalman filter.

In a development of this implementation, the method may further include the steps of, prior to adjusting the actual time of the slave clock, determining the skew of the slave clock according to the updated state equation, and only proceeding to adjust the actual time of the slave clock based on the updated state equation if the skew is less than a predetermined level.

As the slave clock is already frequency synchronized to the master clock via Sync-E, a further advantageous feature of using the Kalman filter model is that it can provide for computation of the frequency offset (α_(n)) as well. The computation of this variable from the Kalman filter can then be used to verify how good the Kalman filtering estimate of the time offset is, that is, when the time offset estimate is good enough to be used for clock alignment by the slave. This is on the basis that if the Kalman filter computes a significant frequency offset when the clocks are known to be frequency synchronised via Sync-E, then it is highly likely that the time offset computation will suffer from similar error and so should not be used to update the local clock in the slave device.

Preferably the step of estimating the variance of the measurement noise calculates a running variance in a recursive manner using the previous estimate of the variance and the calculated measurement noise.

Calculating the variance of a large set of numbers, particularly from real-time data, is not a trivial exercise. For clock data of the type being processed by the method of this aspect, the data set may include millions of numbers. Accordingly, by calculating a running variance in a recursive manner greatly simplifies the process and reduces the processing power required to perform this step. In a particularly preferred embodiment, the estimation of the variance only involves the previous estimate of the variance and the data from the last sample received.

In a development of the method of the present aspect the method may further include the step of, prior to steps of estimating the variance of the measurement noise and adjusting the actual time of the slave clock, screening out data having unusually large measurement noise.

This development of the present aspect provides a packet pre-processing technique (that can be used with the Kalman filtering estimation described above, or any other time synchronization technique) to screen out packets with excessive PDVs which, when used could cause overreaction to measurements.

The method of the present aspect may include any combination of some, all or none of the above described preferred and optional features.

The methods of the above aspect is preferably implemented by a slave device according to the second aspect of this invention or in a system according to the third aspect of this invention, as described below, but need not be.

Further aspects of the present invention include computer for running on computer systems which carry out the methods of the above aspect, including some, all or none of the preferred and optional features of that aspect.

At their broadest, slave devices of the present invention are arranged to synchronise the frequency and time of their local clock with a master clock by using both Synchronous Ethernet communications and additional timing messages.

A second aspect of the present invention provides a slave device connected to a master device having a master clock over an Ethernet network, wherein the slave device includes: a slave clock; an Ethernet adapter; and a processor, wherein: the Ethernet adapter is arranged to receive, at the physical layer, timing signals from the master device via Synchronous Ethernet and the processor synchronises the frequency of the slave clock to the frequency of the master clock using those timing signals; and the Ethernet adapter is arranged to receive, through a layer above the physical layer, first and second timing messages sent from the master device and respective first and second timestamps being the time of sending of the first and second messages according to the master clock, and to record the time of receipt of the first and second messages according to the slave clock; and further wherein the processor is arranged to: calculate the forward variable delay encountered by the first and second messages in travelling from the master device to the slave device over the network by comparing the difference between the interval between the time of sending of said messages and the interval between the time of receipt of said messages; send third and fourth timing messages to the master device for the master device to record the time of receipt of the third and fourth messages according to the master clock; record the time of sending of the third and fourth messages according to the slave clock, receive, from the master device, the recorded time of receipt of the third and fourth messages; calculate the reverse variable delay encountered by the third and fourth messages in travelling from the slave device to the master device over the network by comparing the difference between the interval between the time of sending said messages and the interval between the time of receipt of said messages; calculate, from the forward and reverse variable delays, a measure of the noise (“the measurement noise”) associated with both the receipt of the first and second timing messages and receipt of the third and fourth timing messages; estimate the variance of the measurement noise; and adjust the actual time of the slave clock to synchronize it with the time of the master using the first and third timestamps, the recorded times of receipt of the first and third messages and the estimated variance.

The measure of the noise calculated may be referred to herein as “the measurement noise” but need not be an exact measure of the measurement noise in the system. The estimated variance may be an estimate of the variance of the measure of the noise obtained, which is used as an estimate of the variance of the actual measurement noise, or an estimate of the variance of the actual measurement noise which is produced using the measure of the noise.

The device of this aspect makes use of a new approach for time synchronization where the frequency synchronization features of Synchronous Ethernet (Sync-E) are used by a slave device to extract certain packet transport measurement parameters which can then be used as an input to a time offset computation mechanism used for clock alignment at the slave.

Time synchronization requires both frequency and time offset alignment between a master and slave clock. In the method of the present aspect, frequency synchronization is provided via Sync-E, so the remaining computation is for the slave clock to determine its time offset from the master clock.

The time offset should be fixed (a constant) since slave's frequency is aligned to the master's. However, even for time offset determination, the problem is that timing messages sent end-to-end over the packet network are still exposed to the packet delay variations (PDVs) inherent in the packet network even though frequency synchronization is available through Sync-E. The PDV inherent in packet networks is a primary source of clock noise. The higher the clock noise, the poorer the clock quality rendering the recovered clock sometimes unusable for end system applications when the noise exceeds application defined thresholds.

Preferably the timing messages are messages according to the IEEE 1588 Precision Time Protocol, and more preferably the first and second timing messages are Sync messages and the third and fourth messages are Delay_Req messages.

The IEEE 1588 Precision Time Protocol (PTP) is a method for distributing high precision timing information (both frequency and time) from a master to one or more slaves over a packet network. In a conventional IEEE 1588 PTP arrangement, the slave then uses a clock recovery mechanism to recover the master clock from this timing information. The primary challenge in clock distribution over packet networks is dealing with the variable transit delays experienced by timing packets, that is, the packet delay variations (PDVs). The clock recovery mechanism at a slave will normally have to correct for both frequency and time deviations from the master clock during clock recovery.

Synchronous Ethernet (Sync-E), defined in a suite of ITU-T Recommendations, provides a method for frequency (only) distribution over the physical layer of Ethernet (IEEE 802.3). Since time synchronization requires both frequency and time offset alignment between master and slave clocks, the slave device of the present aspect may combine use of Sync-E and IEEE 1588 PTP where Sync-E provides frequency alignment while IEEE 1588 PTP provides time offset correction.

In this arrangement, Sync-E provides a perfect frequency synchronization signal using Ethernet physical layer frequency transfer, while IEEE 1588 PTP provides time synchronization over the packet network (above the physical layer).

However, the IEEE 1588 PTP messages (or other timing messages) are still exposed to the PDVs inherent in the packet network even though frequency synchronization is available through Sync-E. These PDVs still make time synchronization a challenge even though the frequency synchronization aspect is solved. In preferred implementations of the method of this aspect, time offset alignment with IEEE 1588 PTP with PDV compensation is provided where the synchronized frequency signal available at the slave via Sync-E is used to determine the PDV compensation parameters.

By combining features of Sync-E and IEEE 1588 PTP, the slave device of this aspect may be able to accurately measure the actual PDV of PTP messages arriving at the slave. The slave can then correctly adjust its clock while compensating for the actual delay variations. This may allow the slave to remove the negative effects that these delay variations cause.

In one implementation of the present aspect, the processor adjusts the actual time of the slave clock by using a linear Kalman filter to update a state equation representing the offset and skew of the slave clock compared to the master clock using the estimated variance of the measurement noise following receipt of each timing message in the slave device.

Practical implementation of the Kalman filter typically requires getting a good estimate of the process noise and measurement noise covariances. In the slave device of the present aspect, the frequency synchronization features of Sync-E can be used to extract accurately the PDV parameters from the PTP message flows to the slave. The PDV parameters can then be used to determine the measurement noise covariances for the Kalman filter.

In a development of this implementation, the processor determines the skew of the slave clock according to the updated state equation, and only adjusts the actual time of the slave clock based on the updated state equation if the skew is less than a predetermined level.

As the slave clock is already frequency synchronized to the master clock via Sync-E, a further advantageous feature of using the Kalman filter model is that it can provide for computation of the frequency offset (α_(n)) as well. The computation of this variable from the Kalman filter can then be used to verify how good the Kalman filtering estimate of the time offset is, that is, when the time offset estimate is good enough to be used for clock alignment by the slave. This is on the basis that if the Kalman filter computes a significant frequency offset when the clocks are known to be frequency synchronised via Sync-E, then it is highly likely that the time offset computation will suffer from similar error and so should not be used to update the local clock in the slave device.

Preferably the processor estimates the variance of the measurement noise by calculating a running variance in a recursive manner using the previous estimate of the variance and the calculated measurement noise.

Calculating the variance of a large set of numbers, particularly from real-time data, is not a trivial exercise. For clock data of the type being processed by the method of this aspect, the data set may include millions of numbers. Accordingly, by calculating a running variance in a recursive manner greatly simplifies the process and reduces the processing power required to perform this step. In a particularly preferred embodiment, the estimation of the variance only involves the previous estimate of the variance and the data from the last sample received.

In a development of the present aspect, the processor screens out data having unusually large measurement noise prior to making its estimate of the variance of the measurement noise and adjusting the actual time of the slave clock.

This development of the present aspect provides a packet pre-processing technique (that can be used with the Kalman filtering estimation described above, or any other time synchronization technique) to screen out packets with excessive PDVs which, when used could cause overreaction to measurements.

The slave device of this aspect preferably operates by carrying out the relevant steps of a method according to the above described first aspect.

The slave device of the present aspect may include any combination of some, all or none of the above described preferred and optional features.

At their broadest, systems of the present invention provide for frequency and time synchronisation between a master and slave device by using both Synchronous Ethernet communications and additional timing messages.

A third aspect of the present invention provides a system for synchronising the frequency and time of a slave clock in a slave device to a master clock in a master device, the system including: the master device having the master clock and a first Ethernet adapter; the slave device having the slave clock, a second Ethernet adapter and a processor; and an Ethernet network connecting the master device and the slave device, wherein: the first Ethernet adapter is arranged to send, at the physical layer, timing signals based on said master clock using the Synchronous Ethernet protocol; the second Ethernet adapter is arranged to receive, at the physical layer, said timing signals from the master device and the processor synchronises the frequency of the slave clock to the frequency of the master clock using those timing signals; the master device is arranged to send, through the first Ethernet adapter, using a layer above the physical layer, first and second timing messages and respective first and second timestamps being the time of sending of the first and second messages according to the master clock; the second Ethernet adapter is arranged to receive, through a layer above the physical layer, said first and second timing messages and said first and second timestamps, and to record the time of receipt of the first and second messages according to the slave clock; the slave device is arranged to send, through the second Ethernet adapter, using a layer above the physical later, third and fourth timing messages, and to record the time of sending of the third and fourth messages according to the slave clock; the first Ethernet adapter is arranged to receive, through a layer above the physical layer, said third and fourth timing messages and sent from the slave device and to record the time of receipt of the third and fourth messages according to the master clock and send the recorded time of receipt of the third and fourth messages to the slave device; and further wherein the processor in the slave device is arranged to: calculate the forward variable delay encountered by the first and second messages in travelling from the master device to the slave device over the network by comparing the difference between the interval between the time of sending of said messages and the interval between the time of receipt of said messages; calculate the reverse variable delay encountered by the third and fourth messages in travelling from the slave device to the master device over the network by comparing the difference between the interval between the time of sending said messages and the interval between the time of receipt of said messages; calculate, from the forward and reverse variable delays, a measure of the noise associated with both the receipt of the first and second timing messages and receipt of the third and fourth timing messages; estimate the variance of the measurement noise; and adjust the actual time of the slave clock to synchronize it with the time of the master using the first and third timestamps, the recorded times of receipt of the first and third messages and the estimated variance.

The measure of the noise calculated may be referred to herein as “the measurement noise” but need not be an exact measure of the measurement noise in the system. The estimated variance may be an estimate of the variance of the measure of the noise obtained, which is used as an estimate of the variance of the actual measurement noise, or an estimate of the variance of the actual measurement noise which is produced using the measure of the noise.

The system of this aspect makes use of a new approach for time synchronization where the frequency synchronization features of Synchronous Ethernet (Sync-E) are used by a slave device to extract certain packet transport measurement parameters which can then be used as an input to a time offset computation mechanism used for clock alignment at the slave.

Time synchronization requires both frequency and time offset alignment between a master and slave clock. In the system of the present aspect, frequency synchronization is provided via Sync-E, so the remaining computation is for the slave clock to determine its time offset from the master clock.

The time offset should be fixed (a constant) since slave's frequency is aligned to the master's. However, even for time offset determination, the problem is that timing messages sent end-to-end over the packet network are still exposed to the packet delay variations (PDVs) inherent in the packet network even though frequency synchronization is available through Sync-E. The PDV inherent in packet networks is a primary source of clock noise. The higher the clock noise, the poorer the clock quality rendering the recovered clock sometimes unusable for end system applications when the noise exceeds application defined thresholds.

Preferably the timing messages are messages according to the IEEE 1588 Precision Time Protocol, and more preferably the first and second timing messages are Sync messages and the third and fourth messages are Delay_Req messages.

The IEEE 1588 Precision Time Protocol (PTP) is a method for distributing high precision timing information (both frequency and time) from a master to one or more slaves over a packet network. In a conventional IEEE 1588 PTP arrangement, the slave then uses a clock recovery mechanism to recover the master clock from this timing information. The primary challenge in clock distribution over packet networks is dealing with the variable transit delays experienced by timing packets, that is, the packet delay variations (PDVs). The clock recovery mechanism at a slave will normally have to correct for both frequency and time deviations from the master clock during clock recovery.

Synchronous Ethernet (Sync-E), defined in a suite of ITU-T Recommendations, provides a method for frequency (only) distribution over the physical layer of Ethernet (IEEE 802.3). Since time synchronization requires both frequency and time offset alignment between master and slave clocks, the system of the present aspect may combine use of Sync-E and IEEE 1588 PTP where Sync-E provides frequency alignment while IEEE 1588 PTP provides time offset correction.

In this arrangement, Sync-E provides a perfect frequency synchronization signal using Ethernet physical layer frequency transfer, while IEEE 1588 PTP provides time synchronization over the packet network (above the physical layer).

However, the IEEE 1588 PTP messages (or other timing messages) are still exposed to the PDVs inherent in the packet network even though frequency synchronization is available through Sync-E. These PDVs still make time synchronization a challenge even though the frequency synchronization aspect is solved. In preferred implementations of the method of this aspect, time offset alignment with IEEE 1588 PTP with PDV compensation is provided where the synchronized frequency signal available at the slave via Sync-E is used to determine the PDV compensation parameters.

By combining features of Sync-E and IEEE 1588 PTP, the system of this aspect may be able to accurately measure the actual PDV of PTP messages arriving at the slave. The slave can then correctly adjust its clock while compensating for the actual delay variations. This may allow the slave to remove the negative effects that these delay variations cause.

In one implementation of the present aspect, the processor in the slave device adjusts the actual time of the slave clock by using a linear Kalman filter to update a state equation representing the offset and skew of the slave clock compared to the master clock using the estimated variance of the measurement noise following receipt of each timing message in the slave device.

Practical implementation of the Kalman filter typically requires getting a good estimate of the process noise and measurement noise covariances. In the slave device of the present aspect, the frequency synchronization features of Sync-E can be used to extract accurately the PDV parameters from the PTP message flows to the slave. The PDV parameters can then be used to determine the measurement noise covariances for the Kalman filter.

In a development of this implementation, the processor determines the skew of the slave clock according to the updated state equation, and only adjusts the actual time of the slave clock based on the updated state equation if the skew is less than a predetermined level.

As the slave clock is already frequency synchronized to the master clock via Sync-E, a further advantageous feature of using the Kalman filter model is that it can provide for computation of the frequency offset (α_(n)) as well. The computation of this variable from the Kalman filter can then be used to verify how good the Kalman filtering estimate of the time offset is, that is, when the time offset estimate is good enough to be used for clock alignment by the slave. This is on the basis that if the Kalman filter computes a significant frequency offset when the clocks are known to be frequency synchronised via Sync-E, then it is highly likely that the time offset computation will suffer from similar error and so should not be used to update the local clock in the slave device.

Preferably the processor estimates the variance of the measurement noise by calculating a running variance in a recursive manner using the previous estimate of the variance and the calculated measurement noise.

Calculating the variance of a large set of numbers, particularly from real-time data, is not a trivial exercise. For clock data of the type being processed by the method of this aspect, the data set may include millions of numbers. Accordingly, by calculating a running variance in a recursive manner greatly simplifies the process and reduces the processing power required to perform this step. In a particularly preferred embodiment, the estimation of the variance only involves the previous estimate of the variance and the data from the last sample received.

In a development of the present aspect, the processor screens out data having unusually large measurement noise prior to making its estimate of the variance of the measurement noise and adjusting the actual time of the slave clock.

This development of the present aspect provides a packet pre-processing technique (that can be used with the Kalman filtering estimation described above, or any other time synchronization technique) to screen out packets with excessive PDVs which, when used could cause overreaction to measurements.

The system of this aspect preferably operates by carrying out a method according to the above described first aspect.

The system of the present aspect may include any combination of some, all or none of the above described preferred and optional features.

IEEE 1588 PTP supports both time and frequency synchronization but Sync-E supports high quality frequency synchronization only (at SONET/SDH level quality). The key attributes of the two techniques are given in Table 1 below. Some operators are planning widespread deployment of Sync-E because of its accuracy and SONET/SDH like management capabilities which operators are already very familiar with. The advantages of both PTP and Sync-E have been recognized in the industry and it has been suggested that combining PTP and Sync-E can offer better performance than PTP alone. This allows applications that need time (e.g., TDD wireless technologies) to tap from PTP and those that need frequency (e.g., FDD wireless technologies, circuit emulation services (CES)) to tap from Sync-E. Applications like billing and SLA (service level agreements) can benefit from a network that is supports the time-of-day. The key benefit is both PTP and Sync-E can be used independently or in combination to address the synchronization needs of the network as shown in FIG. 7.

TABLE 1 Sync-E vs. IEEE 1588 PTP Sync-E IEEE 1588 PTP Uses the physical layer of Ethernet Independent of the physical layer Can only distribute frequency - it Can distribute both time-of-day cannot distribute time-of-day and frequency It is not affected by impairments Can be affected by impairments of (PDV, packet loss, out of order) the packet network such as PDV introduced by the higher levels of the network protocols

Development of the Clock Offset and Skew Model

Let θ and α be the clock offset and skew (or normalized frequency deviation), respectively, between a master clock and a slave clock. Below some generalized clock offset and skew equations are set out that will be used in the synchronization solution discussed in more detail later.

At any particular time instant, a snapshot (i.e., instantaneous view) of the relationship between the master (server) clock with timeline S(t) and the slave (client) clock with timeline C(t), can be described by the well-know simple skew clock model depicted in FIG. 8, and described by the equation,

S(t)=(1+α)C(t)+θ,  (1)

where α is a very small quantity in the order of parts-per-million. This snapshot is an instantaneous view of how well the two clocks are (mis)aligned. FIG. 8 explains very well the influence of θ and α on the alignment. The above equation can be extended to account for the case where the master clock and slave clock exchange messages over a communication link with delay.

Assume that a Sync message travelling from a master to a slave experiences a fixed delay d plus a variable (stochastic) delay ε. Similarly, assume that a Delay_Req message sent from the slave to the master experiences a fixed delay of d and a variable delay γ. It is further assumed that the fixed delay components in both directions are equal (symmetric communication paths) but the messages experience variable delays such queuing delays. The master 101 and slave 103 exchange messages using the delay-request delay-response mechanism described in FIG. 6 and also illustrated in FIG. 9. In this system events are logged in discrete time instants n=0, 1, 2, 3, . . . , with each discrete steps is equal to a small interval (sampling) of Δt seconds.

For the nth Sync message (see FIG. 6) which departs with timestamp T_(1,n)εS(t) and arrives with timestamp T_(2,n)εC(t) after having experienced a fixed delay d and a variable delay ε_(n), the simple skew clock model above can be extended to account for the travel time d+ε_(n) to obtain the following expression

(T _(1,n) +d+ε _(n))=(1+α_(n))T _(2,n)+θ_(n)  (2)

The variables θ_(n) and α_(n) are the offset and skew during the nth Sync message exchange.

For the nth Delay_Req message which departs with timestamp T_(3,n) εC(t) and arrives with timestamp T_(4,n)εS(t) after having experienced a fixed delay d and a variable delay γ_(n), the following expression is obtained

T _(4,n) −d−γ _(n)=(1+α_(n))T _(3,n)+θ_(n)  (3)

Adding (2) and (3) gives

T _(1,n) +T _(4,n)+ε_(n)−γ_(n)=(1+α_(n))(T _(2,n) +T _(3,n))+2θ_(n)

(T _(1,n) −T _(2,n))+(T _(4,n) −T _(3,n))=2θ_(n)+α_(n)(T _(2,n) +T _(3,n))+(γ_(n)−ε_(n))  (4)

The expression above can be rewritten as

$\begin{matrix} {\underset{\underset{\gamma_{n}}{}}{\left( {T_{1,n} - T_{2,n}} \right) + \left( {T_{4,n} - T_{3,n}} \right)} = {\underset{\underset{D_{n}X_{n}}{}}{{2\theta_{n}} + {\alpha_{n0}\left( {T_{2,n} + T_{3,n}} \right)}} + \underset{\underset{v_{n}}{}}{\left( {\gamma_{n} - ɛ_{n}} \right)}}} & (5) \end{matrix}$

where n is a nonnegative time index, y_(n)=(T_(1,n)−T_(2,n))+(T_(4,n)−T_(3,n)) is a scalar, D_(n)=[2(T_(2,n)+T_(3,n))] is a vector, X_(n) ^(T)=[θ_(n) α_(n)] is a vector, and ν_(n)=(γ_(n)−ε_(n)) is the measurement noise. The nth sampling interval is considered to be the period in which the nth Sync and nth Delay_Req messages exchanges occur. The clock offset can be expressed from (4) as follows

$\begin{matrix} {\theta_{n} = {{- \left\lbrack \frac{\left( {T_{2,n} - T_{1,n}} \right) - \left( {T_{4,n} - T_{3,n}} \right)}{2} \right\rbrack} - \frac{\alpha_{n}\left( {T_{2,n} + T_{3,n}} \right)}{2} - \frac{\left( {\gamma_{n} - ɛ_{n}} \right)}{2}}} & (6) \end{matrix}$

It can be seen from this equation that for a system with zero skew and no noise, the clock offset simplifies to the classic offset equation

$\begin{matrix} {\theta_{n} = {- \left\lbrack \frac{\left( {T_{2,n} - T_{1,n}} \right) - \left( {T_{4,n} - T_{3,n}} \right)}{2} \right\rbrack}} & (7) \end{matrix}$

The above equation can be derived without the minus sign when the clock model C(t)=(1+α)S(t)+θ is used instead. Either way, the same classic offset equation can be derived.

Discrete-Time Linear Kalman Filter Model

It is assumed that the true process dynamics are described by the following linear transition equation

X _(n) =A _(n) X _(n-1) +B _(n) U _(n) +w _(n),  (8)

where X_(n) is an M-dimensional vector of state variables, A_(n) is an M×M transition matrix, B_(n) is an M×M control-input model matrix which is applied to the M-dimensional control vector U_(n), and w_(n) is an M-dimensional process noise vector which is assumed to be drawn from a zero mean multivariate normal distribution with covariance Q_(n)=E[w_(n)w_(n) ^(T)], w_(n)˜N(0,Q_(n)), and T denotes transpose. In the present estimation problem, the B_(n)U_(n) term is not applicable and so is ignored, resulting in the following state equation

X _(n) =A _(n) X _(n-1) +w _(n)  (9)

The one-step estimate (predicted (a priori) state estimate) is given by

{circumflex over (X)} _(n,n-1) =A _(n) {circumflex over (X)} _(n-1,n-1),  (10)

where, in general {circumflex over (X)}_(n,n-1) is an estimate of X_(n) given data y₁ . . . , y_(n-1) in periods 1 through n−1, where y_(n) is an observed variable in period n. The relations that distinguish the Kalman filter model and associated computational procedures from other linear estimation techniques are the particular model relating X_(n) to y_(n) and the algorithm for computing {circumflex over (X)}_(n,n). The 1×M matrix D_(n) defines the relationship between y_(n) and X_(n) as

y _(n) =D _(n) X _(n)+ν_(n),  (11)

where ν_(n) is the measurement noise which is assumed to be zero mean Gaussian white noise with covariance R_(n)=E[ν_(n)ν_(n) ^(T)], ν_(n)˜N(0,R_(n)). At time n, the vector {circumflex over (X)}_(n,n) is computed by

{circumflex over (X)} _(n,n) ={circumflex over (X)} _(n,n-1) +K _(n) {tilde over (z)} _(n) ={circumflex over (X)} _(n,n-1) +K _(n)(y _(n) −D _(n) {circumflex over (X)} _(n,n-1)).  (12)

This is called the updated (a posteriori) state estimate. The M×1 Kalman gain matrix K_(n) can be calculated recursively by the following equations:

Predicted Error Covariance Matrix:

P _(n,n-1) =A _(n) P _(n-1,n-1) A _(n) ^(T) +Q _(n)  (13)

Innovation (or Residual) Covariance Matrix:

S _(n) =D _(n) P _(n,n-1) D _(n) ^(T) +R _(n)  (14)

Optimal Kalman Gain:

K _(n) =P _(n,n-1) D _(n) ^(T) S _(n) ⁻¹ =P _(n,n-1) D _(n) ^(T) [D _(n) P _(n,n-1) D _(n) ^(T) +R _(n)]⁻¹  (15)

Updated Error Covariance Matrix:

P _(n,n)=(I−K _(n) D _(n))P _(n,n-1)  (16)

It is assumed that E[ν_(n)]=E[w_(n)]=0 for all n, E[w_(n)ν_(i) ^(T)]=0 for all (n,i), E[w_(n)w_(i) ^(T)]=0 for all n≠i, and E[ν_(n)ν_(i) ^(T)]=0 for all n≠i.

From (12), it can be seen that the vector {circumflex over (X)}_(n) which is the estimate of X_(n) is derived as the previous one-step estimation {circumflex over (X)}_(n,n-1), plus a linear combination (weighting) of the differences between the measurement y_(n) and the previous estimate of these measurements D{circumflex over (X)}_(n,n-1). The weights assigned to the difference terms are appropriate components of the gain matrix K_(n). It is also important to note that {circumflex over (X)}_(n,n) depends on {circumflex over (X)}_(n,n-1), y_(n), K_(n), and D, but not explicitly on y₁ . . . , y_(n-1). This recursive nature of the Kalman filter eliminates the need for storage of historical data.

It can be seen from (15) that K_(n) has terms which are directly proportional to the elements of the covariance matrix Q and inversely proportional to the elements of R. That is, K_(n) is, in a sense, proportional to the variability of the true process dynamics and inversely proportional to the measurement variability. Thus, it is the “Q/R” relationship which defines the responsiveness of the filter (via the K matrix) to estimated errors in state (y_(n)−D_(n){circumflex over (X)}_(n,n-1)). As the elements of K_(n) decrease (by decreasing Q or increasing R), the estimates become more stable. That is, if there is more confidence in the unbiasedness of the process model or less confidence in the observations, the filter will be designed to respond more slowly to apparent deviations from the predicted trend line. As K_(n) increases (by increasing Q or decreasing R), estimates detect and respond to data that deviate from the assumed linear model, represented by the expected value of (8).

Thus, in summary the estimation procedure has the following steps:

-   -   1. When n=0, the filter is initialized by a user-supplied         estimate {circumflex over (X)}_(0,0) of the initial state vector         X₀ and P_(0,0)     -   2. Using these estimates, a one-step estimate is produced using         {circumflex over (X)}_(n,n-1)=A_(n){circumflex over         (X)}_(n-1,n-1)     -   3. The Kalman gain matrix K_(n) and the matrix P_(n,n) are         calculated     -   4. When a new observation is received, {circumflex over         (X)}_(n,n)={circumflex over         (X)}_(n,n-1)+K_(n)(y_(n)−D_(n){circumflex over (X)}_(n,n-1)) is         used to obtain an estimate of the present state vector     -   5. Then using this new estimate {circumflex over (X)}_(n,n), go         to Step 2 and repeat process.

The algorithm continues to process new observations and produce estimates in this manner.

There are several points to be noted:

-   -   1. If the matrices D, R, A, and Q are independent of n, the gain         matrix K_(n) and the matrices P and S are independent of the         observations. Thus, these matrices can be pre-calculated for use         in the algorithm.     -   2. No past observations must be stored since all historical         information is contained in the estimate {circumflex over         (X)}_(n,n) or equivalently via {circumflex over         (X)}_(n,n-1)=A_(n){circumflex over (X)}_(n-1,n-1), the one-step         ahead estimate.     -   3. Only the one-step ahead estimate must be saved to be used in         the next step's process.

When the assumptions listed above about the process and measurement noises are valid and the models accurately describe the true process dynamics and the measurement system, the Kalman filter produces unbiased estimates of {circumflex over (X)}_(n,n-1), that is, E[{circumflex over (X)}_(n,n-1)]=X_(n,n-1). In addition, the estimates {circumflex over (X)}_(n,n-1) have minimum variance in the class of all unbiased estimators.

New Technique for Characterizing the Network PDV and System Measurement (Observation) Noise Using the Sync-E Clock

If the all the Kalman filter model matrices (A, D, Q, and R) are known, the algorithm is completely specified; moreover, the desirable statistical properties of the filter are assured. However, in general, the true model is not known precisely. Usually, it is the case that the processes being modeled suggest a reasonable choice of matrices A and D. Also, typically, the data characteristics can be analyzed so that Q and R are known analytically or empirically. In the present clock synchronization problem, the D matrix is known from (5) as D_(n)=[2 (T_(2,n)+T_(3,n))].

With the reference frequency available at both master and slave through the Sync-E communication link, the slave has only to correct for time offset in the face of PDV. This is still a challenging task even though the frequency alignment aspect of the problem has been solved by using Sync-E between master and slave.

From (5) and assuming the measurement noise satisfies the properties given above, the covariance of the measurement noise ν_(n)=(γ_(n)−ε_(n)) can be expressed as

$\begin{matrix} {R_{n} = {{E\left\lbrack {v_{n}v_{i}^{T}} \right\rbrack} = \left\{ \begin{matrix} {\sigma_{v,n}^{2} = \sigma_{{({\gamma - ɛ})},n}^{2}} & {n = i} \\ 0 & {otherwise} \end{matrix} \right.}} & (17) \end{matrix}$

which is simply the variance of the measurement error or equivalently the variance of the difference between the Delay_Req message PDVs, γ_(n) and the Sync message PDVs, ε_(n). Unlike Q which is a matrix, R is not, it is just a scalar quantity.

Below a method according to an embodiment of the present invention is set out for characterizing ε_(n) and γ_(n) using the PTP message flows and the Sync-E clock available at both master and slave.

When Sync messages are transmitted over the packet network, they will arrive at the slave with variable delay, the PDV, ε. The Sync messages are generated using the high quality Sync-E clock available at the master. This clock is typically traceable to a Primary Reference Clock (PRC), the high quality reference clock of the network. Any arbitrary Sync message can be taken as an initialization (or anchor) message as shown in FIG. 10. The timestamp difference between the nth and (n−1)th generated Sync message at the master is defined as ΔT_(f,n)=T_(1,n)−T_(1,n-1). The timestamp difference between the nth and (n−1)th Sync message arrivals as measured by the Sync-E clock available at slave is defined as ΔR_(f,n)=T_(2,n)−T_(2,n-1). Since Sync-E clock is the same as at the master, the timestamp difference measured by the slave includes the PDV, ε_(n), experienced between the two arrivals, that is, ε_(n)=ΔR_(f,n)−ΔT_(f,n).

The above can be applied to Delay_Req messages transmitted from the slave to the master. The timestamp difference between the nth and (n−1)th generated Delay_Req message at the slave can be defined as ΔT_(r,n)=T_(3,n)−T_(3,n-1). The timestamp difference between the nth and (n−1)th Delay_Req message arrivals as measured by the Sync-E clock available at master is defined as ΔR_(r,n)=T_(4,n)−T_(4,n-1). The timestamp difference measured by the master includes the PDV, γ_(n), experienced between the two arrivals, that is, γ_(n)=ΔR_(r,n)−ΔT_(r,n). The calculation of the ε_(n) and γ_(n) from the PTP messages is illustrated in FIG. 11.

The measurement noise can then be obtained from the forward PDV ε_(n) and the reverse PDV γ_(n) as ν_(n)=(γ_(n)−ε_(n)). The covariance of ν_(n), which simplifies to the variance σ_(ν,n) ², has to be determined to be used in the Kalman filter estimation of the clock offset. Calculating the variance of a large set of numbers, particularly from real-time data, is not trivial. It becomes more challenging when data is generated in the order of millions of numbers and even beyond. Using the naive method for calculating variance would therefore not be efficient. In the present method, the calculation of the running variance using a recursion relationship involving only the last sample is used as follows.

The following definitions are used: ν _(n) as the running mean of ν_(n) up to discrete time n, με(0,1) as a weighting parameter, and V_(n) as the running variance of ν_(n) at time n. The running mean can be calculated using the well-known exponentially weighted moving average (EWMA) as

ν _(n)=(1−μ) ν _(n-1)+μν_(n)  (18)

In general, the variance of a finite data of size n is given by

$\begin{matrix} \begin{matrix} {V_{n} = {\frac{1}{n}{\sum\limits_{i = 1}^{n}\left( {v_{i} - {\overset{\_}{v}}_{n}} \right)^{2}}}} \\ {= {\frac{1}{n}\left\lbrack {\left( {v_{n} - {\overset{\_}{v}}_{n}} \right)^{2} = {\sum\limits_{i = 1}^{n - 1}\left( {v_{i} - {\overset{\_}{v}}_{n}} \right)^{2}}} \right\rbrack}} \\ {= {\frac{1}{n}\left\lbrack {\left( {v_{n} - {\overset{\_}{v}}_{n}} \right)^{2} = {\left( {n + 1} \right)\frac{1}{\left( {n - 1} \right)}{\sum\limits_{i = 1}^{n - 1}\left( {v_{i} - {\overset{\_}{v}}_{n}} \right)^{2}}}} \right\rbrack}} \end{matrix} & (19) \end{matrix}$

Now by setting

${V_{n - 1} = {\frac{1}{\left( {n - 1} \right)}{\sum\limits_{i = 1}^{n - 1}\left( {v_{i} - {\overset{\_}{v}}_{n}} \right)^{2}}}},$

in the above equation, we get

$\begin{matrix} {V_{n} = {{\frac{1}{n}\left\lbrack {\left( {v_{n} - {\overset{\_}{v}}_{n}} \right)^{2} + {\left( {n - 1} \right)V_{n - 1}}} \right\rbrack} = {{\frac{1}{n}\left( {v_{n} - {\overset{\_}{v}}_{n}} \right)^{20}} + {\left( {1 - \frac{1}{n}} \right)V_{n - 1}}}}} & (20) \end{matrix}$

Since n, is a positive number, 1/n will have to be a constant between 0 (if n is infinite) and 1 (if n=1). Substituting β for 1/n, we have

V _(n)=β(ν_(n)− ν _(n))²+(1β)V _(n-1),  (21)

where the weighting factor βε(0,1). The weights μ and β can be set equal if that works well for the application. At each iteration of the Kalman filter estimation, the running variance V_(n) is determined and then ρ_(ν,n) ²=V_(n) is set for use in the Kalman filter.

FIG. 9 shows the distribution and processing of timing information in a slave device 3 according to an embodiment of the present invention. The PTP protocol engines 11, 21 in the master and slave 20 implement the delay request-response mechanism and exchange PTP messages according to FIG. 6 for the two-step clock case. Both protocol engines are driven by Sync-E clocks traceable to the frequency reference 1. The protocol engine 21 at the slave 20 outputs a set of timestamps {T1, T2, T3, T4} after each execution of the delay request-response mechanism. With a current and previous sets of timestamps {T1, T2, T3, T4}, a processor in the slave 20 extracts the relevant noise statistics (23, 24) for the Kalman filter estimation of the clock offset (25) as described above. Using the estimated clock offset and the Sync-E driven local clock 22, a processor in the slave 20 estimates the server (master) time (26).

Characterizing the Clock (Process) Noise

In this section the clock (process) model parameters A and w_(n) are derived. The A matrix can be obtained as follows. The clock skew between two points T_(1,n) and T_(1,n-1) can be estimated given two clock offsets θ_(n) and θ_(n-1) as

$\begin{matrix} {\alpha_{n - 1} = {\frac{\theta_{n} - \theta_{n - 1}}{T_{1,n} - T_{1,{n - 1}}}.}} & (22) \end{matrix}$

The process dynamics for the clock while accounting for process noise can then be expressed as

θ_(n)=θ_(n-1)+α_(n-1)(T _(1,n) −T _(1,n-1))+w _(θ,n)

α_(n)=α_(n-1) +w _(α,n)  (23)

where w_(n) ^(T)=[w_(θ,n), w_(α,n)] is the process noise vector which is assumed to be drawn from a zero mean normal distribution with covariance Q_(n)=E[w_(n)w_(n) ^(T)]. The system can be described by the following two-state dynamic model

$\begin{matrix} {{X_{n} = {\begin{bmatrix} \theta_{n} \\ \alpha_{n} \end{bmatrix} = {{{\begin{bmatrix} 1 & \left( {T_{1,n} - T_{1,{n - 1}}} \right) \\ 0 & 1 \end{bmatrix}\begin{bmatrix} \theta_{n - 1} \\ \alpha_{n - 1} \end{bmatrix}} + \begin{bmatrix} w_{\theta,n} \\ w_{\alpha,n} \end{bmatrix}} = {{AX}_{n - 1} + w_{n}}}}},} & (24) \end{matrix}$

where A_(n) is the known 2-by-2 state transition matrix. If the time between Sync messages is fixed, then, ΔT_(n)=(T_(1,n)−T_(1,n-1))=Δt is a constant term, and we have

$\begin{matrix} {A = {\begin{bmatrix} 1 & \left( {T_{1,n} - T_{1,{n - 1}}} \right) \\ 0 & 1 \end{bmatrix} = {\begin{bmatrix} 1 & {\Delta \; t} \\ 0 & 1 \end{bmatrix}.}}} & (25) \end{matrix}$

The process noise covariance Q is related to the oscillator of local PLL at the slave (FIG. 7). For the Telecom clock synchronization problem, the local oscillator is typically a Temperature Compensated Crystal Oscillator (TCXO) or Oven Controlled Crystal Oscillator (OCXO).

Typically, the clock deviation of a clock is modeled by the following random differential equation

$\begin{matrix} {{\frac{{X(t)}}{t} = {{{FX}(t)} + {\xi (t)}}},} & (26) \end{matrix}$

where F is an L×L real matrix and X(t) is the state vector

$\begin{matrix} {{X(t)} = {\begin{bmatrix} {x_{1}(t)} \\ \vdots \\ {x_{L}(t)} \end{bmatrix}.}} & (27) \end{matrix}$

The quantity x₁(t) represents the time deviation (or time offset) and x₂(t) is part of the clock frequency deviation. In the model, ξ(t) is the vector of statistically independent zero-mean white Gaussian noises,

$\begin{matrix} {{{\xi (t)} = \begin{bmatrix} {\xi_{1}(t)} \\ \vdots \\ {\xi_{L}(t)} \end{bmatrix}},} & (28) \end{matrix}$

with each element of the vector having autocorrelation

E[ξ _(k)(t ₁)ξ_(k)(t ₂)]=g _(k)δ(t ₁ −t ₂),  (29)

for k=1, 2, . . . , L and where δ(t) is the Dirac delta function. The autocorrelation matrix of ξ(t) is thus

E[ξ(t ₁)ξ^(T)(t ₂)]=Gε(t ₁ −t ₂),  (30)

and where G is the diagonal matrix

$\begin{matrix} {G = {\begin{bmatrix} g_{1} & \ldots & 0 \\ \vdots & \ddots & \vdots \\ 0 & \ldots & g_{L} \end{bmatrix}.}} & (31) \end{matrix}$

The random differential (26) can be written in the integral form

X(t)=Φ(t−t ₀)X(t ₀)+∫_(t) ₀ ^(t)Φ(t−τ)ξ(τ)dτ,  (32)

where t₀ is the initial time and Φ(t) is the transition matrix

Φ(t)=e ^(Ft).  (33)

The time axis can be discretised with a sampling time Δt so that t takes the values t_(n)=nΔt, where n=0, 1, 2, . . . . The integral for of (32) in discrete time becomes

X(t _(n))=Φ(Δt)X(t _(n-1))+w(t _(n-1)),  (34)

where

w(t _(n-1))=∫_(t) _(n-1) ^(t) ^(n) Φ(t−τ)ξ(τ)dτ,  (35)

is a Gaussian random variable whose mean is zero and whose covariance matrix is

$\begin{matrix} \begin{matrix} {Q_{n - 1} = {Q\left( t_{n - 1} \right)}} \\ {= {E\left\lbrack {\left( {{w\left( t_{n - 1} \right)} - {E\left\lbrack {w\left( t_{n - 1} \right)} \right\rbrack}} \right)\left( {{w\left( t_{n - 1} \right)} - {E\left\lbrack {w\left( t_{n - 1} \right)} \right\rbrack}} \right)^{T}} \right\rbrack}} \\ {= {\int_{t_{n - 1}}^{t_{n}}{\int_{t_{n - 1}}^{t_{n}}{{\Phi \left( {t_{n} - \tau} \right)}{E\left\lbrack {{\xi (\tau)}{\xi^{T}(s)}} \right\rbrack}{\Phi^{T}\left( {t_{n} - s} \right)}\ {s}\ {\tau}}}}} \\ {= {\int_{t_{n - 1}}^{t_{n}}{\int_{t_{n - 1}}^{t_{n}}{{\Phi \left( {t_{n} - \tau} \right)}G\; {\delta \left( {s - \tau} \right)}\ {\Phi^{T}\left( {t_{n} - s} \right)}{s}\ {\tau}}}}} \\ {= {\int_{t_{n - 1}}^{t_{n}}{{\Phi \left( {t_{n} - \tau} \right)}G\; {\Phi^{T}\left( {t_{n} - \tau} \right)}\ {\tau}}}} \end{matrix} & (36) \end{matrix}$

For L=2, we note that the matrix A in our clock process model (25) is

$\begin{matrix} {{{A \approx {\Phi \left( {\Delta \; t} \right)}} = {{I + {F\; \Delta \; t}} = \begin{bmatrix} 1 & {\Delta \; t} \\ 0 & 1 \end{bmatrix}}},{F = \begin{bmatrix} 0 & 1 \\ 0 & 0 \end{bmatrix}},} & (37) \end{matrix}$

which an approximation of a power series expansion of (33), which is

$\begin{matrix} {{{\Phi (t)} = {^{Ft} = {I + {F(t)} + {F^{2}\frac{t^{2}}{2}} + {F^{3}\frac{t^{2}}{6}} + \ldots + {F^{q - 1}\frac{t^{q - 1}}{\left( {q - 1} \right)!}}}}},} & (38) \end{matrix}$

where I is the identity matrix. With L=2, and F defined in (37), the time deviation of the clock is modeled by a linear deterministic trend plus a Brownian motion and an integrated Brownian motion. From (26) we have

$\begin{matrix} {{\frac{{x_{1}(t)}}{t} = {{x_{2}(t)} + \eta_{1} + {g_{1}{\xi_{1}(t)}}}}{\frac{{x_{2}(t)}}{t} = {\eta_{2} + {g_{2}{\xi_{2}(t)}}}}} & (39) \end{matrix}$

where we account for deterministic trends by choosing the deterministic nonzero initial conditions, represented by the constants, η₁ and η₂.

$\begin{matrix} {{X(0)} = {\begin{bmatrix} {x_{1}(0)} \\ {x_{2}(0)} \end{bmatrix} = \begin{bmatrix} \eta_{1} \\ \eta_{2} \end{bmatrix}}} & (40) \end{matrix}$

The constants can be interpreted as the deterministic phenomena driving the clock errors.

Substituting A and G (in (31) for L=2), we have

$\begin{matrix} \begin{matrix} {Q_{n - 1} = {Q\left( {\Delta \; t} \right)}} \\ {= {\int_{0}^{\Delta \; t}{{\Phi \left( {{\Delta \; t} - \tau} \right)}G\; {\Phi^{T}\left( {{\Delta \; t} - \tau} \right)}{\tau}}}} \\ {= {\int_{0}^{\Delta \; t}{{{\begin{bmatrix} 1 & \left( {{\Delta \; t} - \tau} \right) \\ 0 & 1 \end{bmatrix}\begin{bmatrix} g_{1} & 0 \\ 0 & g_{2} \end{bmatrix}}\begin{bmatrix} 1 & 0 \\ \left( {{\Delta \; t} - \tau} \right) & 1 \end{bmatrix}}\ {\tau}}}} \\ {= \begin{bmatrix} \left( {{g_{1}\Delta \; t} + \frac{g_{2}\Delta \; t^{3}}{3}} \right) & \frac{g_{2}\Delta \; t^{2}}{2} \\ \frac{g_{2}\Delta \; t^{2}}{2} & {g_{2}\Delta \; t} \end{bmatrix}} \end{matrix} & (41) \end{matrix}$

The parameters g₁ and g₂ play the following role in the clock noise process. The parameter g₁ drives the white frequency noise, corresponding to a random walk on the time deviation (offset), while g₂ drives the frequency random walk, corresponding to an integrated random walk on the time deviation. The parameters g₁ and g₂ can also be denoted as the spectral amplitudes of f₁ and f₂, respectively, for the process model in (24). For white noise, the spectral amplitudes g₁ and g₂ can be expressed in terms of well know Allan variance parameters. The Allan variance of the frequency deviation f(t) of a clock can be computed as [1][2][3]

$\begin{matrix} {{{\sigma_{f}^{2}(\tau)} = {\frac{1}{2}{\langle\left( {{\overset{\_}{f}\left( {t + \tau} \right)} - {\overset{\_}{f}(t)}} \right)^{2}\rangle}}},} & (42) \end{matrix}$

where τ>0 is the observation interval, the operator



denotes the time average

$\begin{matrix} {{{\langle\left( {{\overset{\_}{f}\left( {t + \tau} \right)} - {\overset{\_}{f}(t)}} \right)^{2}\rangle} = {\lim\limits_{T\rightarrow\infty}{\frac{1}{T}{\int_{{- T}/2}^{T/2}{\left( {{\overset{\_}{f}\left( {t + \tau} \right)} - {\overset{\_}{f}(t)}} \right)^{2}\ {t}}}}}},} & (43) \end{matrix}$

and f(t) is the average frequency defined by

$\begin{matrix} {{\overset{\_}{f}(t)} - {\frac{1}{\tau}{\int_{t - \tau}^{t}{{f(s)}\ {s}}}}} & (44) \end{matrix}$

The average frequency value and the increment ( f(t+τ)− f(t)) are therefore functions of τ, in which case the Allan variance σ_(f) ²(τ) also depends on τ. We note that the Allan variance of the frequency deviation f(t) of a cesium clock can be computed as [1][2][3]

$\begin{matrix} {{{\sigma_{f}^{2}(\tau)} = {\frac{g_{1}}{\tau} + {\frac{g_{2}}{3}\tau}}},} & (45) \end{matrix}$

where term containing g₁ is the Allan variance of the white frequency noise, while the term containing g₂ corresponds to the random walk frequency noise.

The Allan variance is finite for all clock noises and has the power-law behavior [4]

σ_(f) ²(τ)=K _(σ)τ^(μ),  (46)

where K_(σ) is a constant. Taking the logarithm of this:

log σ_(f) ²=log K _(σ)+μ log τ,  (47)

which indicates that the Allan variance of power-law noise is a straight line in bilogarithmic coordinates. Every clock has a different signature in the Allan variance domain, which can be used to classify the different types of clocks and noises. By estimating the slopes of the Allan variance it is possible to identify the types of noise in the clock signal, and to understand possible anomalies in the stability performances.

The Allan variance is the accepted standard for the stability analysis of oscillators and is commonly used by clock manufacturers to specify the performance of their products. By taking measurable quantities such as time differences, a clock can be compared against a reference clock (with superior performance), and then the estimated deviations and noises due to the lower performance clock found. The parameters g₁ and g₂ can be determined from the Allan variance of the client clock from which the process covariance matrix Q can be calculated.

Time Offset Accuracy Improvement with PTP Message Pre-Processing—Outlier Screening

In this section a technique is set out which screens out unusual data or outliers in the measurement process described by (5). “Unusual data” here is defined as that representing an unusually large (outlier) measurement error (ν_(n)). The importance of these types of errors is that the estimation process tends to adjust the previous estimate {circumflex over (X)}_(n,n-1) in the direction of the new measurement y_(n), with the movement being proportional to the estimated error (y_(n)−D_(n){circumflex over (X)}_(n,n-1)). Therefore, without additional logic of an adaptive nature, outliers would cause overreaction to measurements. A schematic layout of the interaction of this pre-processing with the other timing synchronization features within a PTP slave device is shown in FIG. 12. The PTP slave device in FIG. 12 is essentially identical to that previously described and shown in FIG. 9, except that the PTP slave device 20 has a PTP Message Pre-processing function 27 which operates on all PTP Timestamps extracted from and applied to the PTP messages to screen out those which are considered unreliable.

Measurements are considered to be carried out over a measurement period (window) Δt where a number of Sync and Delay_Req messages are exchanged between the master and slave. Note that Δt here is not the time between Sync message transmissions (as described in (23) to (25)) but rather a window consisting of a number of Sync and Delay_Req transmissions. One of the following two methods can be used to select the “best” (that is, optimal or “less contaminated”) Sync message and Delay_Req message among all messages exchanged in a measurement period.

1. Select a Sync Message and Delay_Req Message Belonging to the Same Protocol Exchange Pair:

-   -   Here, a representative Sync message plus its associated         Delay_Req message belonging to the same protocol exchange pair         are selected and their timestamps used in the estimation process         of both PDV and offset elimination (as shown in FIG. 12) in         order to arrive at the estimate of server time. The running         variance V_(n) described above (see Equation (21)) for the         variable ν_(n)=(γ_(n)−ε_(n)) can be used in the selection         process. The goal here is to select the message pair that         satisfies a screening criterion, that is, the message pair whose         PDV has the minimum variance among all other message pairs in         the measurement period.

2. Select a Sync Message and Delay_Req Message not Necessarily Belonging to the Same Protocol Exchange Pair:

-   -   Out of the exchanges in a window, one representative Sync         message and one Delay_Req message is selected and their         timestamps used in the estimation process (as shown in FIG. 12).         The goal here is to select the messages that satisfy a screening         criterion, that is, messages that have PDVs that carry the         minimum variance among all other messages. The PDVs experienced         by messages (ε_(n) and γ_(n)) are determined according to the         procedure explained above in relation to FIG. 10 and FIG. 11. A         running variance of the PDV ε_(n) is calculated for each Sync         message in a sampling interval.

ν _(ε,n)=(1−μ) ν _(ε,n-1)+με_(n)  (48)

V _(ε,n)=β(ε_(n)− ν _(ε,n))²+(1−β)V _(ε,n-1),  (49)

Similarly a running variance of the PDV γ_(n) is determined for each Delay_Req message.

ν _(γ,n)=(1−μ) ν _(γ,n-1)+μγ_(n)  (50)

V _(γ,n)=β(γ_(n)− ν _(γ,n))²+(1−β)V _(γ,n-1),  (51)

To minimize storage of parameters, the computation of the variance is done sequentially in real-time as messages arrive in such a way that at the end of a sampling interval, only the Sync and Delay_Req messages that resulted in the minimum variance among all messages are kept. There is no need to compute running variances for all messages and then select the minimum from them. It is assumed that the client oscillator is stable enough within the sampling interval as is the case with quartz oscillators. Fortunately, a quartz oscillator (a TCXO and OCXO) is an excellent choice when short-term stability is needed. At the end of a sampling internal, a Sync message and Delay_Req message with associated timestamps are used in the estimation process.

Time Offset Computation

It is noted that the Kalman filter model (24) gives both time offset θ_(n) and frequency offset α_(n) estimates. However, in the present clock synchronization problem using Sync-E and PTP, frequency synchronization is already achieved (given) by Sync-E, and as a result, we expect α_(n) to converge to a small number, in the order of the accuracies of a Stratum 1 clock (10⁻¹¹). The results for α_(n) is a bonus that we get from our Kalman filter model—since the frequency offset α_(n) is part of the state vector X, we get a frequency estimate α_(n) along with the time offset estimate.

This quantity α_(n) can be used to verify how good the Kalman filtering estimate of the time offset is, that is, α_(n) must converge to a very small number of this order. The converged time offset value is then used to adjust the client clock to obtain and estimate of the server clock as illustrated in FIG. 13.

The systems and methods of the above embodiments may be implemented in a computer system (in particular in computer hardware or in computer software) in addition to the structural components and user interactions described.

The term “computer system” includes the hardware, software and data storage devices for embodying a system or carrying out a method according to the above described embodiments. For example, a computer system may comprise a central processing unit (CPU), input means, output means and data storage. Preferably the computer system has a monitor to provide a visual output display (for example in the design of the business process). The data storage may comprise RAM, disk drives or other computer readable media. The computer system may include a plurality of computing devices connected by a network and able to communicate with each other over that network.

The methods of the above embodiments may be provided as computer programs or as computer program products or computer readable media carrying a computer program which is arranged, when run on a computer, to perform the method(s) described above.

The term “computer readable media” includes, without limitation, any medium or media which can be read and accessed directly by a computer or computer system. The media can include, but are not limited to, magnetic storage media such as floppy discs, hard disc storage media and magnetic tape; optical storage media such as optical discs or CD-ROMs; electrical storage media such as memory, including RAM, ROM and flash memory; and hybrids and combinations of the above such as magnetic/optical storage media.

While the invention has been described in conjunction with the exemplary embodiments described above, many equivalent modifications and variations will be apparent to those skilled in the art when given this disclosure. Accordingly, the exemplary embodiments of the invention set forth above are considered to be illustrative and not limiting. Various changes to the described embodiments may be made without departing from the spirit and scope of the invention.

In particular, although the methods of the above embodiments have been described as being implemented on the systems of the embodiments described, the methods and systems of the present invention need not be implemented in conjunction with each other, but can be implemented on alternative systems or using alternative methods respectively.

REFERENCES

-   [1]. A. J. Van Dierendonck and R. G. Brown, “Relationship between     Allan variances and Kalman filter parameters,” in Proc. 16th Annual.     Precise Time and Time Interval (PTTI) Meeting, Greenbelt, Md., 1984,     pp. 273-292. -   [2]. C. Zucca and P. Tavella, “The clock model and its relationship     with the Allan and related variances,” IEEE Trans. Ultrason.,     Ferroelect., Freq. Contr., vol. 52, no. 2, pp. 289-296, 2005. -   [3]. W. Chaffee, “Observability, ensemble averaging and GPS time,”     IEEE Trans. Aerosp. Electron. Syst., vol. 28, no. 1, pp. 224-240,     1992. -   [4]. D. W. Allan, “Statistics of atomic frequency standards,” Proc.     IEEE, vol. 54, no. 2, pp. 221-230, 1966.     All references referred to above are hereby incorporated by     reference. 

1. A method of synchronising the frequency and time of a slave clock in a slave device to a master clock in a master device, wherein the master device and the slave device are connected by an Ethernet network, the method including the steps of: synchronising the frequency of the slave clock to that of the master clock using Synchronous Ethernet communications between the master device and the slave device; receiving in the slave device first and second timing messages sent from the master device and respective first and second timestamps being the time of sending of the first and second messages according to the master clock, and recording the time of receipt of the first and second messages according to the slave clock; calculating the forward variable delay encountered by the first and second messages in travelling from the master device to the slave device over the network by comparing the difference between the interval between the time of sending of said messages and the interval between the time of receipt of said messages; receiving in the master device third and fourth timing messages sent from the slave device and recording the time of receipt of the third and fourth messages according to the master clock; recording in the slave device the time of sending of the third and fourth messages according to the slave clock; sending the recorded time of receipt of the third and fourth messages to the slave device; calculating the reverse variable delay encountered by the third and fourth messages in travelling from the slave device to the master device over the network by comparing the difference between the interval between the time of sending said messages and the interval between the time of receipt of said messages; calculating, from the forward and reverse variable delays, a measure of the noise associated with both the receipt of the first and second timing messages and receipt of third and fourth timing messages; estimating the variance of the measurement noise; and adjusting the actual time of the slave clock to synchronize it with the time of the master using the first and third timestamps, the recorded times of receipt of the first and third messages and the estimated variance.
 2. A method according to claim 1 wherein the timing messages are messages according to the IEEE 1588 Precision Time Protocol, and the first and second timing messages are Sync messages and the third and fourth messages are Delay_Req messages.
 3. A method according to claim 1 wherein the step of adjusting the actual time of the slave clock uses a linear Kalman filter to update a state equation representing the offset and skew of the slave clock compared to the master clock using the estimated variance of the measurement noise following receipt of each timing message in the slave device.
 4. A method according to claim 3, further including the steps of, prior to adjusting the actual time of the slave clock, determining the skew of the slave clock according to the updated state equation, and only proceeding to adjust the actual time of the slave clock based on the updated state equation if the skew is less than a predetermined level.
 5. A method according to claim 1 wherein the step of estimating the variance of the measurement noise calculates a running variance in a recursive manner using the previous estimate of the variance and the calculated measurement noise.
 6. A method according to claim 1 further including the step of, prior to steps of estimating the variance of the measurement noise and adjusting the actual time of the slave clock, screening out data having unusually large measurement noise.
 7. A slave device connected to a master device having a master clock over an Ethernet network, wherein the slave device includes: a slave clock; an Ethernet adapter; and a processor, wherein: the Ethernet adapter is arranged to receive, at the physical layer, timing signals from the master device via Synchronous Ethernet and the processor synchronises the frequency of the slave clock to the frequency of the master clock using those timing signals; and the Ethernet adapter is arranged to receive, through a layer above the physical layer, first and second timing messages sent from the master device and respective first and second timestamps being the time of sending of the first and second messages according to the master clock, and to record the time of receipt of the first and second messages according to the slave clock; and further wherein the processor is arranged to: calculate the forward variable delay encountered by the first and second messages in travelling from the master device to the slave device over the network by comparing the difference between the interval between the time of sending of said messages and the interval between the time of receipt of said messages; send third and fourth timing messages to the master device for the master device to record the time of receipt of the third and fourth messages according to the master clock; record the time of sending of the third and fourth messages according to the slave clock; receive, from the master device, the recorded time of receipt of the third and fourth messages; calculate the reverse variable delay encountered by the third and fourth messages in travelling from the slave device to the master device over the network by comparing the difference between the interval between the time of sending said messages and the interval between the time of receipt of said messages; calculate, from the forward and reverse variable delays, a measure of the noise associated with both the receipt of the first and second timing messages and receipt of third and fourth timing messages; estimate the variance of the measurement noise; and adjust the actual time of the slave clock to synchronize it with the time of the master using the first and third timestamps, the recorded times of receipt of the first and third messages and the estimated variance.
 8. A slave device according to claim 7 wherein the timing messages are messages according to the IEEE 1588 Precision Time Protocol, and the first and second timing messages are Sync messages and the third and fourth messages are Delay_Req messages.
 9. A slave device according to claim 7 wherein the processor adjusts the actual time of the slave clock by using a linear Kalman filter to update a state equation representing the offset and skew of the slave clock compared to the master clock using the estimated variance of the measurement noise following receipt of each timing message in the slave device.
 10. A slave device according to claim 9 wherein the processor determines the skew of the slave clock according to the updated state equation, and only adjusts the actual time of the slave clock based on the updated state equation if the skew is less than a predetermined level.
 11. A slave device according to claim 7 wherein the processor estimates the variance of the measurement noise by calculating a running variance in a recursive manner using the previous estimate of the variance and the calculated measurement noise.
 12. A slave device according to claim 7 wherein the processor screens out data having unusually large measurement noise prior to making its estimate of the variance of the measurement noise and adjusting the actual time of the slave clock.
 13. A system for synchronising the frequency and time of a slave clock in a slave device to a master clock in a master device, the system including: the master device having the master clock and a first Ethernet adapter; the slave device having the slave clock, a second Ethernet adapter and a processor; and an Ethernet network connecting the master device and the slave device, wherein: the first Ethernet adapter is arranged to send, at the physical layer, timing signals based on said master clock using the Synchronous Ethernet protocol; the second Ethernet adapter is arranged to receive, at the physical layer, said timing signals from the master device and the processor synchronises the frequency of the slave clock to the frequency of the master clock using those timing signals; the master device is arranged to send, through the first Ethernet adapter, using a layer above the physical layer, first and second timing messages and respective first and second timestamps being the time of sending of the first and second messages according to the master clock; the second Ethernet adapter is arranged to receive, through a layer above the physical layer, said first and second timing messages and said first and second timestamps, and to record the time of receipt of the first and second messages according to the slave clock; the slave device is arranged to send, through the second Ethernet adapter, using a layer above the physical later, third and fourth timing messages, and to record the time of sending of the third and fourth messages according to the slave clock; the first Ethernet adapter is arranged to receive, through a layer above the physical layer, said third and fourth timing messages and sent from the slave device and to record the time of receipt of the third and fourth messages according to the master clock and send the recorded time of receipt of the third and fourth messages to the slave device; and further wherein the processor in the slave device is arranged to: calculate the forward variable delay encountered by the first and second messages in travelling from the master device to the slave device over the network by comparing the difference between the interval between the time of sending of said messages and the interval between the time of receipt of said messages; calculate the reverse variable delay encountered by the third and fourth messages in travelling from the slave device to the master device over the network by comparing the difference between the interval between the time of sending said messages and the interval between the time of receipt of said messages; calculate, from the forward and reverse variable delays, a measure of the noise associated with both the receipt of the first and second timing messages and receipt of third and fourth timing messages; estimate the variance of the measurement noise; and adjust the actual time of the slave clock to synchronize it with the time of the master using the first and third timestamps, the recorded times of receipt of the first and third messages and the estimated variance.
 14. A system according to claim 13 wherein the timing messages are messages according to the IEEE 1588 Precision Time Protocol, and the first and second timing messages are Sync messages and the third and fourth messages are Delay_Req messages.
 15. A system according to claim 13 wherein the processor adjusts the actual time of the slave clock by using a linear Kalman filter to update a state equation representing the offset and skew of the slave clock compared to the master clock using the estimated variance of the measurement noise following receipt of each timing message in the slave device.
 16. A system according to claim 15 wherein the processor determines the skew of the slave clock according to the updated state equation, and only adjusts the actual time of the slave clock based on the updated state equation if the skew is less than a predetermined level.
 17. A system according to claim 13 wherein the processor estimates the variance of the measurement noise by calculating a running variance in a recursive manner using the previous estimate of the variance and the calculated measurement noise.
 18. A system according to claim 13 wherein the processor screens out data having unusually large measurement noise prior to making its estimate of the variance of the measurement noise and adjusting the actual time of the slave clock. 